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Abstract 

We compute the damping rates of one-particle excitations in a cold ultra- 
relativistic plasma to leading order in the coupling constant e for three types 
of interaction: Yukawa coupling to a massless scalar boson, QED and QCD. 
Damping rates of charged particles in QED and QCD are of order e 3 /U, while 
damping rates of other particles are of order e 4 /i or e 4 /xlog(l/e). We find 
that the damping rate of an electron or of a quark is constant far from the 
Fermi surface, and decreases linearly with the excitation energy close to the 
Fermi surface. This unusual behavior is attributed to the long-range magnetic 
interactions. 
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I. INTRODUCTION 



The quasiparticle concept is a powerful tool for studying the dynamical properties of ul- 
trarelativistic plasmas; it has been widely used in the recent litterature Weakly excited 
states of a plasma can be described as superpositions of quasiparticle states, which behave, 
at least in a first approximation, as free, non-interacting, particles. These elementary exci- 
tations undergo damping through their mutual coupling, which gives their energy spectrum 
a finite width. The quasiparticle concept holds only if this energy spread, or equivalently 
the damping rate, is negligible compared to the mean energy of their wave packet. 

In this paper, we compute damping rates in plasmas at high density and zero temper- 
ature, T = 0. While the thermodynamic properties of ultrarelativistic degenerate plasmas 
are well known J2|, little work has been devoted so far to their dynamical properties HIJ. 
Most works on damping rates have been focusing on plasmas at high temperature with 
vanishing chemical potential, fi = 0. However, plasmas with finite chemical potential are 
also relevant in view of phenomeno logical applications: first, degenerate quark systems (with 
/i<T) might exist in the cores of neutron stars || ; second, the state of the hadronic matter 
(possibly a quark-gluon plasma) temporarily formed in an ultrarelativistic nucleus-nucleus 
collision corresponds to values of \i and T of the same order of magnitude at the presently 
available energies ||. 

From a theoretical point of view, the zero temperature case owes its interest to the basic 
difference, compared to the high temperature limit, in the infrared divergences which plague 
perturbative calculations in gauge theories. Naive perturbation theory leads to two levels 
of infrared divergences J7|,|8| . Those appearing in the calculation of transport coefficients 



and collisional energy loss [T(| are cured by taking into account screening effects at the 
one-loop level: propagators and vertices must be corrected at long wavelengths, according 
to the resummation scheme developed by Braaten and Pisarski [pT[| . The momentum and 
energy dependence of these medium effects is the same for hot (yU = 0) and cold (T = 0) 
relativistic plasmas, and is characterized by a scale of order eT for a hot plasma, and e\i for 
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a cold plasma. At high temperature, the lack of static screening of the transverse part of 
the interaction is responsible for a logarithmic divergence in the perturbative calculation of 
the fermion damping rate |T2|-|l4||. This problem has been recently solved in hot QED ||15 



By contrast, this second level of divergence does not appear at zero temperature, because 
of Pauli blocking. This allows us to obtain finite, most often analytical expressions for the 
damping rates T(p) of all one-particle excitations, as a function of their momentum p. 

Three types of interactions are discussed and compared. (1) Electromagnetic interaction 
(QED): this corresponds to a degenerate electron gas. (2) Non-abelian SU(N C ) interaction 
(QCD): for N c = 3, this corresponds to degenerate quark matter. (3) In addition, we consider 
a toy model where the fermion field i/j(x) is coupled to a massless scalar field (f)(x), with 
the Yukawa interaction Cy = eijjifxf). These theories share several features. The resummed 
fermion propagator is the same for the three theories. The resummed boson propagator is 
the same for the two gauge theories, QED and QCD, while screening corrections reduce to a 
medium-induced mass term |16| for the Yukawa interaction. The coupling constant, noted 
e for all theories, will be assumed to be much smaller than unity, e < 1, to ensure that a 
perturbative expansion is reliable. It has been shown that perturbation theory can be used 
to study the properties of the QCD phase if the temperature and/or the chemical potential 
is much larger than Aqcd @- 

We are considering degenerate Fermi systems in their ground state. Single particle 
excitations are obtained by adding or removing one particle from the ground state, and 
damping of these excitations results from collisions with the fermions of the Fermi sea. Two 
types of excitations compose the quasiparticle spectrum. At large ("hard") momenta p, 
p ~ /i ^> e/i, the elementary modes correspond to single particles slightly perturbed by the 
medium. Their damping processes are the object of this paper. For hard fermion or hole 
excitations, the only collisional process to leading order in e is elastic (M0ller) scattering 
// — ► //. [We use the generic names bosons (b), fermions (/) and antifermions (/).] For an 
antifermion, two processes may contribute, namely elastic (Bhabha) scattering // — > ff and 
pair annihilation // — > bb. For a boson, the only process is elastic (Compton) scattering 



bf — > 6/. On the other hand, long wavelength ("soft") modes, of momentum p of order 
e/i, correspond to collective oscillations in the medium. Their damping processes are more 
complex: to leading order in e, bremsstrahlung contributes: // — > //& for a soft hole, or 
/// — > 6/ for a soft antifermion. A discussion of these calculations will be presented in a 
forthcoming publication JT7J. 

Depending on the interaction and on the type of excitation considered, the dominant 
contribution to the damping rate may come either from processes with large scattering 
angles (9 ~ 1), from processes with small scattering angles (9 ~ e), or from both. Large 
scattering angles correspond to hard exchanged quanta, for which medium effects are small: 
then, the amplitudes are the same as in the vacuum. On the other hand, if a hard particle 
of momentum p ~ \x is deflected by an angle of order e, the exchanged quantum has a 
soft momentum of order e/x, for which medium effects become important. In section ||, 
we estimate the orders of magnitude of both hard and soft contributions to the damping 
rates. We show how they can be calculated, in a kinetic approach, as phase space integrals of 
elementary scattering processes. The equivalence with field theoretical methods is recalled in 
Appendix B. Section |TJ presents the calculations of damping rates for the various excitations 
considered, while section [TV] is devoted to a discussion of the results along with a comparison 
between damping rates in cold and hot plasmas. 

II. GENERAL ANALYSIS 

A. Orders of magnitude 

The damping rate T, the number of collisions per unit time, is of order V ~ anv where 
er is the scattering cross section, n the density of scatterers and v the relative velocity. In 
this paper, we assume that all particles are massless, which implies v = 1 (H = c = 1). 
We consider 2 — > 2 elementary processes in which an incoming particle of four-momentum 
P = (p,p) (fermion above the Fermi level, antifermion or boson) is added to the system 
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and scatters on a particle of the Fermi sea with four-momentum K = (k,k), k < fi. We 
denote by P' = (p',p') and K' = (k',k f ) the four-momenta of outgoing particles. If the 
excitation under study is a hole in the Fermi sea, incoming and outgoing particles must 
be interchanged. The differential cross section for unpolarized particles can generally be 
written as 

"i = S^W, (i) 

where \M\ 2 denotes the scattering matrix element squared, averaged over the helicity states 
of the incoming particle with momentum p, and summed over the helicity states of the 
other particles. (Note that we do a sum, rather than an average, over the 2 helicity states 
of the scatterer. This unusual convention will turn out to be convenient in the following 
sections.) A factor 1/2 is included in \M\ 2 if the outgoing particles are identical. The 
variables s = (P + K) 2 , t = (P — P') 2 and u = (P — K') 2 = —s — t are the usual Mandelstam 
variables. In the center of mass frame, t is related to the momentum transfer q by t = — q 2 
and to the scattering angle 9 by t = — s sin 2 (#/2). In order to obtain the total cross section, 
one must integrate eq. ([TJ) over t between t = — s and t = 0. 

For a given scattering process in the plasma, the squared matrix element \M\ 2 is deduced 
from the usual Feynman rules with appropriate corrections taking into account medium 
polarization effects. We distinguish three cases, depending on the behaviour of \M\ 2 at 
small t: 

1) The tree matrix element \M\ 2 is finite at small momentum transfers: one example is 
fermion-fermion scattering in a Yukawa theory, for which the tree matrix element is 

|M|2 = 3e 4 . (2) 

The total cross section is therefore finite and of order a ~ e 4 //x 2 , since s is of order fi 2 . The 
density of particles per unit volume in the Fermi sea is of order n ~ fi 3 , which gives T ~ e \x. 

2) \M\ 2 is proportional to 1/t 2 at small t: the differential cross section grows like 1/8 4 at 
small scattering angle 6, which makes the total cross section diverge. This occurs in collisions 
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between charged particles (Rutherford divergence). For electron-electron scattering in QED 
(M0ller scattering), 



so that the total cross section diverges linearly in the infrared: da ~ e^dt/t 2 . However, 
in a plasma, due to the screening of the electric charge, the interaction potential decreases 
exponentially for distances larger than the Debye screening length r D ~ l/(e/i). This causes 
a saturation of the differential cross section for momentum transfers q < r^ 1 ~ e/i, i.e. 
for scattering angles 6 < e, or equivalently for t < t min = e 2 /i 2 . The total cross section is 
therefore of order a ~ e 4 /t min ~ e 2 //i 2 . Note that this argument is not completely valid 
for the transverse part of the interaction, which is screened dynamically only (i.e. at finite 
energy transfer). However, as we shall show in detail later, dynamical screening is sufficient 
at zero temperature to saturate the damping rate at small momentum transfer. Now, for 
small momentum transfers of order e/i, the density of scatterers is no longer n ~ /i 3 : because 
of Pauli blocking, the outgoing electron must have an energy larger than the Fermi energy 
/i, and only the electrons just below the Fermi surface can contribute, within an interval of 
order e/i. Their density is n ~ e/i 3 , and one finally obtains r ~ an ~ e 3 /i. 

3) |M| 2 is proportional to 1/t at small t. This happens in processes where the exchanged 
particle is a massless fermion, such as pair annihilation and Compton scattering. For in- 
stance, in the Yukawa theory, the square tree matrix element for pair annihilation into two 
massless Yukawa bosons (// — > bb) is 



and the total cross section diverges logarithmically in the infrared, rather than linearly in the 
previous case. This divergence is also cured by collective effects, which become important 
when the momentum transfer is of order e/i. Note that these medium effects are taken into 
account in the fermion propagator, instead of the boson propagator in the previous case. 
Following the same reasoning, the total cross section is thus of order a = e 4 ln(l/e)//x 2 . 




(3) 




(4) 
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For an annihilation process, there is no Pauli blocking in the final state, and the density of 
scatterers is of order n ~ /i 3 , hence the damping rate is of order V ~ e 4 ln(l/e)/x. In the 
case of Compton scattering, this last remark applies only for a boson whose energy is larger 
than the Fermi energy, so that it can turn into a fermion with almost the same momentum. 
Otherwise, Pauli blocking alone inhibits the infrared divergence and the damping rate is of 
order e 4 /i (section |Tl [Pi) . 



These considerations on orders of magnitude allow us to simplify the damping rate cal- 
culation. In case 1), kinematics is dominated by hard momentum transfer (the "hard sector 
contribution"): the phase space for soft momentum transfer (the "soft sector contribution") 
is smaller by a factor e. On the other hand, soft momentum transfer dominates in case 2); 
the differential cross section integrated over momentum transfer of order e/i is larger by a 
factor 1/e 2 than the contribution from hard momentum transfer, while the phase space is 
smaller by only a factor e. Finally, in case 3), both hard and soft sector contributions are 
of equal magnitude. 

We separate the hard and soft scales by introducing an arbitrary IR (resp. UV) momen- 
tum cutoff in the hard (resp. soft) contribution. The cutoff g* is chosen at an intermediate 
scale, e/j C g, <C fx. The hard and soft sector contributions to the damping rate are both 
computed in the framework of kinetic theory. For q > (hard contribution), collective 
effects are negligible since q 3> e/i, and the damping rate can be calculated from the vacuum 
scattering amplitudes, integrated over phase space with an appropriate choice of variables, 
described in section |II B| . For q < q* (soft contribution), medium effects must be taken into 
account, but the kinematics is simplified by the fact that the exchanged particle is much 
softer than the external ones. The damping rate is most simply calculated as the emis- 
sion probability of a soft, off-shell quantum, as developed in section [11 Q In case 3), both 



approximations are compatible in the region of the cutoff, so that the total damping rate, 
which is the sum of the two contributions, does not depend on the cutoff g* [ p~8f| . 

Note that processes involving soft external particles, either incoming or outgoing, are 
subleading because the associated phase space is small. Thus, the external legs in the 



processes considered in this paper will always be hard. The situation is different, of course, 
when studying the damping of soft excitations |I7|| , where at least the incoming particle is 
soft. 



B. Hard sector contribution 

The hard sector contribution is calculated as the total transition rate, integrated over 
the available phase space: 

r h (p) = ^J dr p/ dT k dr kt jM^(27c) A 5\P + K - P' - K'). (5) 

Here, dr k = d 3 k/2k(2ir) 3 is the Lorentz invariant phase space volume. When all the particles 
involved in the scattering process are hard, including the exchanged one, medium effects can 
be ignored to leading order, and the matrix element is computed with the usual Feynman 
rules. If the outgoing particles are electrons, the phase space is limited by the Pauli exclusion 
principle: their energy must be larger than the Fermi energy \i. 

Alternatively, we could have chosen to compute the damping rate in a field theoretical 
approach, as the imaginary part of the self-energy. Generally, the leading hard sector con- 
tribution corresponds to the imaginary part of two-loop energy diagrams, as illustrated in 
Fig. 1 in the case of electron-electron scattering. This equivalence is proven explicitly in 
Appendix B. 

The phase space integration in eq. (§) is easily carried out if an appropriate choice of 
kinematic variables is made: there are 9 integration variables, and 4 constraints, hence 
5 degrees of freedom. Since the transition probability is invariant under a simultaneous 
rotation of p', k and k' about the direction of p, by integrating over this angle only 4 
degrees of freedom remain. At this point, it is convenient to introduce the 4- momentum 
Q = (u, q) transferred by the incoming particle to the plasma: 

Q = (co,q) = P-P' = K'~K (6) 



S 



We choose 3 variables as k, uj and q = \q\. Note that k and uj fix the energy of the 
incoming and outgoing particles, while q fixes the angle between k and k' . The fourth 
degree of freedom corresponds to the angle between the plane spanned by (p,p r ) and the 
plane spanned by (k, k') or, equivalently, to the azimuthal angle of k around q. With 
these variables, eq. @ becomes 

T h {p) = J ^Jdkdujdq(\M?) (7) 

where the brackets ( ) denote an average over the azimuthal angle 0. 

The integration limits on q are easily derived from the definitions in eq. @: 

\uj\ < q < min(A; / + k,p' + p) = min(2A; + uj, 2p — uj). (8) 

Since the scattered fermion belongs to the Fermi sea, we have the constraint k < fi. Depend- 
ing on whether or not the outgoing particles are fermions, additional constraints on k and 
uj may result from the Pauli blocking conditions: k! = k + uj > \x and/or p' = p — uj > \x. Fi- 
nally, the momentum transfer is kept hard by imposing q > q* with e/i < g» <C //. Together 
with eq. (g), these conditions completely specify the integration domain in eq. (0). 

C. Soft sector contribution 

When the momentum transfer of a scattering process is of order e/z, one must correct 
the bare interaction for medium (screening) effects, whose contribution in the propagator 
is of the same order as the bare propagator itself. This is achieved by replacing the bare 
propagator by the resummed propagator given in Appendix A. The damping rate is then 
evaluated as a collision integral similar to eq. (|5|), however with a screened interaction, as 
illustrated in Fig. 2 (left) for electron-electron scattering. 

The soft contribution to the damping rate can in fact be cast into a much simpler form, as 
the transition rate of a process where the the incoming hard particle, with four-momentum 
P = (p,p), emits a virtual (spacelike) soft particle with four-momentum Q = (uj,q) and 



9 



scatters into a hard particle with four-momentum P' = (p',p f ) (Fig. 2, right). The fact that 
the soft particle is spacelike means that the actual physical process is, in the example of 
Fig. 2, electron-electron scattering, not "electron decay". We show in Appendix B that the 
two approaches (Fig. 2, left, and Fig. 2, right) are equivalent and amount to evaluate the 
imaginary part of a (resummed) one-loop self energy diagram (Fig. 2, middle). 

The damping rate V(p), for the process in Fig. 2 (right), can be evaluated from Fermi's 
Golden Rule, as in eq. (g|), with an important modification: the soft energy u and momentum 
q — \q\ are no longer related by a dispersion relation. Instead, a whole range of values are 
allowed according to a spectral distribution whose actual form is derived from the screening 
corrections. The damping rate takes then the following form: 

Tsip) = ^Jj^y piQ) 27n5((p " Q)2) (9) 

In this expression, p{Q) denotes the spectral function of the soft particle. The spectral 
function of the outgoing hard particle is simply 2n5((P — Q) 2 )- \M\ 2 denotes the matrix 
element of the transition on the right of Fig. 2, squared and summed over final spins, and 
averaged over the spin states of the incoming particle. 

In the limit where q is much softer than p, (P — Q) 2 = —2p(u — qcos8), where 9 is 
the angle between q and p. Thus the condition that the outgoing hard particle is on its 
mass shell, (P — Q) 2 = reduces to cos 9 — u/q. Integrating over 9, the previous equation 
becomes 

= TTtW f* ^ f df>P(u,q)W- (10) 

lblT z p z JO J-q 

Note that u> > (resp. uj < 0) corresponds to the emission (resp. absorption) of a soft 
particle. 

The phase space is further restricted by the Fermi-Dirac and Bose-Einstein distributions, 
eqs. (|B5|). If the outgoing hard particle is a fermion, it must be above the Fermi level, which 
implies u < p — \i. If the soft particle is a boson, the only possibility is u > because 
there is no boson initially present at T = 0. If the soft particle is a fermion, both signs are 
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possible: 00 < corresponds to the absorption of a soft fermion, and u > to the emission 
of a soft antifermion. Note that the phase space for the soft virtual particle is limited by 
the same statistical constraints as if it were a real, on-shell, particle. 

It turns out that the matrix element squared \M\ 2 always takes a very simple form, as 
we shall see in the next section. 



III. RESULTS 



We now compute explicitly the damping rates of one-particle excitations to leading order 
in e for the three theories (Yukawa, QED and QCD). We first study the fermionic excitations, 
in sections |III A| to |III C[ then the bosonic excitations, in sections [1.1 ID] and |IIIE| . For each 
type of excitation, we show the self-energy diagrams whose imaginary part corresponds to 
the elementary process under consideration. 



A. Fermion in Yukawa's theory 

We consider an incoming fermion of momentum p above the Fermi level: p > fi. In the 
Yukawa theory, the fermion-fermion scattering matrix element squared, corresponding to 
the Feynman diagrams depicted in Fig. 3, is a constant given by eq. (0). Collisions are thus 
dominated by hard momentum transfers, as discussed in section [II A| . The rate T is given 
by eq. (]?]), where Fermi statistics imposes the conditions 

0</jL — k<u<p — fx. (11) 

Using eq. (|TT]), the limits on q given by eqs. (|[) reduce to u < q < 2k + u>. The integration 



for p > 2ji, 
for fi < p < 2ji. 



is straightforward and leads to the expressions 

e 4 /1 2 (3^0-4^) 



r / x _ I 1287T 3 p 2 no x 

e 4 G0-/1) 2 (ifi-p) U - ) 



128tt 3 p 



2 
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This damping rate corresponds to the imaginary part of the two-loop diagrams displayed in 
Fig. 4. The third two-loop diagram, the rainbow diagram (see Fig. 7) does not contribute at 
zero temperature: its imaginary part corresponds to Compton scattering or pair annihilation. 

We now consider a hole state (of momentum p, p < jj). The scattering process is now 
described in two steps: first an initial vacancy in the Fermi sea is filled by a fermion of 
momentum p'; next the energy difference p' — p is transferred to a fermion of momentum 
k', which is then extracted out of the Fermi sea. Hence, Fermi statistics imposes for this 
process the conditions < k',p' < fi < k. We define uj and q as in eq. (H). Thus eqs. (^) 



and (H) are still valid. However, eq. (jTT|) is now replaced by 



p — ji < u < jj, — k < 0. (13) 

With these conditions, eq. (^) reduces to — uj < q < 2p — uj. The phase space integration of 
eq. (|7|) gives then 

r(p) = {fi ~ p)2 . (14) 

yF! 128tt 3 p v ; 

Near the Fermi surface, the damping rate of fermions and holes vanishes quadratically with 
the excitation energy \p — y\. We shall come back to this in section 0. 

Note that the damping rate diverges for small p. Extrapolating the above formula to the 
soft domain p ~ e/i (where our calculation does not apply), one guesses that Y is of order 
e 3 /i for a soft excitation, instead of e 4 /x for a hard excitation. A correct calculation shows 



that it is indeed the case [17] 



B. Antifermion in Yukawa's theory 

Two collision processes contribute to leading order : Bhabha scattering (Fig. 5 a) and 
pair annihilation (Fig. 5 b). 

The first one gives no difficulty Its matrix element is 



|M| 2 = 2e 4 , (15) 
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and the phase space integration goes along the same lines as in section [III A| , except for 



the fact that no restriction applies on the final state energy p'. Equation ( [TT]) is therefore 
replaced by 

0<fi-k<u. (16) 

The contribution of Bhabha scattering to the damping rate is then obtained by integrating 
eq. (0) using eqs. (H) and ([IB]) : 

e 4 

r ^) = r 9 e"V ( 17 ) 
l 192^^ (3p-2/i) forp> ^ 
For pair annihilation, the tree matrix element is given by eq. (f|), and we need to include 

both hard and soft momentum transfers, according to the discussion following this equation. 

The hard sector contribution is given by eq. (0). To use this equation, we must average the 

matrix element over 0, the azimuthal angle of k with respect to q. We first note that 

exchanging t and u in eq. @ amounts to exchanging the two outgoing photons, thus the u/t 



and t/u terms give identical contributions and \M\ 2 can be replaced by e A (u/t — 1). From the 
definition of t and eq. @, t = Q 2 = uj 2 — q 2 . The variable u is given by u = 2(p ■ k' — pk'). 
Decomposing p and k' into longitudinal and transverse components with respect to q, 
and averaging over 0, one gets (p ■ k') = (p ■ q){k' • q)/q 2 . From eq. @, one obtains 
p ■ q = up + (q 2 — uj 2 )/2 and k' ■ q = uk + (q 2 + uj 2 )/2. The average over gives therefore 
(|M| 2 ) = e\(u/t) - 1), with 

(j) = ^- 2 [(2k + uj)(2p-uj)-q 2 }. (18) 

Notice that the collinear divergence (at uj ~ ±g) canceled out of the ratio. 

The limits on q are given by eq. (§), and the only additional restriction from Fermi 
statistics is k < /i. The integral is logarithmically divergent: 



1287r d p V } ' 



where g* is an IR cutoff for the q— integral. 
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The soft sector contribution corresponds to pair annihilation "at low angles" , in the sense 
that the outgoing photons have momenta very close to those of the incoming electron and 
positron (Fig. 6, left). This can be viewed as a process where a hard antifermion "turns" 
into a hard boson (see Fig. 6, right) by absorbing a soft fermion (u < 0) or emitting a soft 
antifermion (u > 0). The spectral function of a soft fermion receives contributions from 
two channels which are labeled by + and — in appendix A. We denote by M + and M_ the 
corresponding matrix elements 

M± = eu{p,\)u{±q,\'). (20) 

Only states with opposite chiralities have non-vanishing matrix elements, which implies 
A' = —A. Using eq. (|A7|) and the property that 



cos (0/2), denoting the angle 



between p and q, one easily obtains the result 

£|M ± | 2 = 2e 2 p(g^), (21) 



A' 



where we have used the relation cos# = u/q. From eqs. ( ID ) and Q2"I|), the soft contribution 
to the damping rate becomes 

e 2 fi* fi 

Ts ^ > = ^ r plo qdq J_ q du [( q ~ u }P+( UJ > q ) + ( q + UJ ")P-( UJ > q ^> ( 22 ) 
or, by using the relationship p + (ui,q) = p-(—u,q), 

e 2 [i* [i 

F 2s (p) = J Q qdq J did (q - u) p + {q,u). (23) 

To integrate over u, we first show the following sum rule: 

du (q - u)p + (q, uj) = 0. (24) 

Since p + is the discontinuity of the Green function G + on the real axis, the contour integral 
in eq. (p4]) is simply the integral of (q — z)G+(q, z), with zona contour going parallel and 
right above the real axis from — oo to +oo, and coming back right below it from +oo to — oo. 
Deforming the contour into a circle of infinite radius, on which G+(q, z) reduces to the free 
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propagator [2q(z — q)}" 1 , one easily sees that the integral vanishes. With the help of the 
sum rule of eq. (|24|), one can then express the contribution from the cut piece (— q < uo < q) 
of the density p + in eq. (p3|) in terms of the pole piece (\u\ > q) of the density p+. The pole 
contribution to the spectral density is p + = 2n5(G^ 1 ). Using eqs. ( |A8| ) and (|A9|) , we get 



r 2s (p) 



47r 2 pJo qdq l\ |> duJ ( u ~ ^^fai^ ~ q) 



(25) 



Since the mass operator H + (q,ui), given by eq. ( |A 1 (J| ) , depends only on u/q, it is convenient 
to change variables from q, u to x = u/q and y = 2q{uj — q), which gives 



t 2s (p) 



ydxdy 



5 y 



(26) 



IGnp J \x — 1| 

This can be readily integrated over y and then over x. The integration limits are derived 
from the dispersion relation eq. (|A11|) . For q ^> rrif, one branch is at x ~ 1 + mj / q 2 and 
the other is at x ~ — 1. Thus the integral on x extends from — oo to —1 and from 1 + mj/q 2 
to +oo. One gets 



2™2 



r 2s (p) 



em 



f 



inp 



m 



f. 



(27) 



Replacing m/ by its value given in Appendix A and adding eqs. (flTf ) and (^), the cutoff 
cancels out. The contribution of pair annihilation to the damping rate is thus 



T 2 (p)=T 2h (p) + T 2s (p) 

( 2 VP\ 



eV 
128tt 3 j9 



log 



\m 2 f J 



(28) 



The total damping rate to order e 4 is the sum of Ti and T 2 given by Eqs. ( |TTD and 
(P5|). It corresponds to the imaginary part of the self-energy diagrams displayed in Fig. 7. 
The two-loop diagrams give the hard contribution, while the (resummed) one-loop diagram 
gives the soft contribution. 



C. Fermion and antifermion in QED and QCD 



We turn now to studying the fermion lifetime in relativistic QED and QCD plasmas. 
The damping process is fermion-fermion scattering and is dominated by soft momentum 
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transfers, as explained in section [II B| . In QED, we compute this process as the emission of 
an off-shell soft photon by a hard electron (Fig. 2). The matrix element of the transition is 
M = eu(p', A')7 M w(p, A)e M (q), where e^q) is the polarization of the photon. Since q <C p, 
one can replace p' = p — q by p in the matrix element. Then, the first factor in M is 
simply the electric current associated with the incoming electron, which reduces to J M = 
eu(p, A')7 M w(p, A) = 2eP^5\ t y ) the Kronecker symbol reflecting the fact that helicity is 
conserved in the process. We thus obtain M = 2e P ■ e. Now, in the Coulomb gauge, the 
spectral function of the soft photon receives a contribution from longitudinal and transverse 
modes, which are denoted by pz,(g, uj) and pr(q,uj) respectively (see Appendix A). The 



longitudinal polarization vector gives a matrix element \Ml\ 2 = Ae 2 p 2 . The two transverse 
polarization vectors e^(q, A) satisfy J2\=i,2 6 t(Qi ^) e T(<?> ^) = — lilj/l 2 which gives, upon 



using uj = qcosO, \M T \ 2 = 4e 2 p 2 (l — u 2 /q 2 ). The soft contribution to the damping rate is 
given by eq. ( |10D and the condition < uj < p — \i: 

( „>2\ 1 

(29) 



e 2 fl* /■min(g,p-/i) 

r a (p) = qdq du 



o 



p L (u, q)+(l- j p T (uj, q) 



Ait 2 Jo 

The restriction uj > is due to the conditions k < p and k' > p, which together imply 
uj = k' — k > (see Fig. 2, left). The damping of a hole is calculated in the same manner 
and leads to the same expression eq. (|29|) , with p — p replaced by p — p. 

Since the spectral functions pi and px fall off rapidly for q > ep, one can safely take 
the cutoff g* to infinity. Introducing then the dimensionless quantities y = q/qo in eq. (f29|). 
e = \p — p\/qD, and x = uj/q, one obtains 

r ( £ ) = -4-2 J dx J y 2d v [ r L(yi x ) + (! - x 2 ) r T(y, x)\ , (30) 

with r LjT = q 2 D p LjT . 

In QCD, the damping rate of a quark or a quark hole has a similar expression because 
the color current of a quark (which enters the matrix element) has the same structure as 
its electric current. The result is then given by eq. (|30| ) with a multiplicative color factor 
Cf = (N 2 — 1)/ (2iV c ): N 2 — 1 for the number of soft gluons, 1/N C for the average over the 
quark colors and 1/2 for the trace of SU(3) generator products. 
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The integrals in eq. ( j30|) can be evaluated numerically. Simple approximate results can 
be derived in two limits: 

1. Far from the Fermi surface, e 1. As the spectral densities ri and ry vanish rapidly 
for y > 1, we can extend the upper bound of the y— integral to infinity. One finds 

r(p) = 0M7e 2 C f q D , (31) 

or 

f 0.018e 3 /i for QED, 

r(p) = r- , (32) 

[ 0M7^N f e 3 /i for QCD with JV) flavors, N c = 3. 

2. Very close to the Fermi surface, e<l. The integration over y covers an appreciable 
range of values for small values of x ~ e only, so that one can safely extend the integral 
over x to infinity. In the static limit, x = uo/q <C 1, the longitudinal and transverse spectral 
functions behave very differently. Static electric fields are screened at distances larger than 
q^ 1 . The longitudinal polarization function, given by eq. ( |A4| ), reduces to a constant, II L ~ 
—qjji and the corresponding spectral function tl is 

ri = PTTF (33) 

Integrating over x first, one finds a contribution to the damping rate of order e 2 



Tl = C f ^£. r dy f £/y dx ^l— 2 + 0(, 3 ) ~ C f ^e 2 + 0(e 3 ). (34) 



On the contrary, the transverse polarization function is purely imaginary for i C 1, because 
a static magnetic field is not screened. In this limit (i = w/g < 1), eq. ( |A5|) gives II t = 
—inqfyUj/Aq 3 , and the spectral function tt is 

TT 2(y 4 + 7r 2 (x/4) 2 )' 1 ' 

Even though there is no static screening, the term proportional to x 2 in the denominator 
induces a deviation from the Rutherford 1/g 4 (i.e. 1/y 4 ) behavior, which is referred to as 
dynamical screening. The main contribution to the damping rate comes from values of y 
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and x (x ~ e) such that y 4 and x 2 are of the same order, i.e. such that q ~ (g^a;) 1 / 3 . Upon 
introducing the variable y = y'y/x/2 and integrating over x first, one finds 

'% £ jf*'p^ +0 ^~ 'w e + ^- (36) 

The damping rate is thus dominated by the transverse contribution and gives 

f^-b-^l forQED, 
r(p) = ^ (37) 
I — b - /u| for QCD with JV) flavors, Y c = 3. 

97T 

For the damping rate of a positron or an antiquark, Rutherford scattering has a matrix 
element varying as {u/t) 2 and dominates over pair-annihilation and the s— channel contri- 
bution to Bhabha scattering, varying respectively as (u/t) and (u 2 /s 2 ). The soft sector 
contribution gives therefore the leading order term. Since the electric and color currents are 
the same as for electrons or quarks, up to a sign, the damping rate is also given by eq. fl30| ) . 
without the restriction t < e/x coming from Pauli blocking. The damping rate is then given 
by the result of eq. for all momenta. 

D. Boson in Yukawa's theory and QED 

An incident beam (of given central energy) of scalar bosons or photons will undergo 
a spectral broadening due to the elastic scattering of its quanta with electrons (Compton 



scattering). The matrix element squared, |M| 2 , corresponding to the diagrams depicted in 
Fig. 8, is related to that of pair annihilation by crossing symmetry. For a Yukawa interaction, 
one deduces from eq. (|J) 



|M| 2 = 2e 4 (---- + 2), (38) 

S XL 



whereas for a hard photon we have 



|M| 2 = 4e 4 (----). (39) 
s u , 



The computation of the damping rate is similar to the calculation done in section [D FBI. In 
particular, both hard and soft momentum transfers may contribute. 



To compute the hard contribution, we define the variables uo = p + k = p' + k! and 



q = p + k = p' + k', so that s 



— q 2 . Taking k', to and q as integration variables, 



the contribution to the damping rate is given by eq. (|7|), with k replaced by k'. The matrix 



element |M| 2 must be averaged over the azimuthal angle of k around q. Following the same 
method as for pair annihilation, we obtain an equation similar to eq. (|PS|): 

1 



2q 2 



{2k' - u){2p - u) -q' 



(40) 



The limits on q are max(|2p— u\, \2k'—u\) < q < cj, and the constraints from Fermi statistics 
are u — p < \i <kl '. 

The inverse term (—s/u) is integrated using the variables u = k — p' = k' — p, q = 
k — k' = k' — p and following the same steps as for the direct term. 

Adding up the direct and inverse terms, we find for the hard sector contribution to the 
damping rate of a Yukawa boson 



T{p < fx) 
T(p > Li) 



For a photon (QED) we obtain 

r(p < n) 

T(p > Li) 



4 

e p 

1927T 3 

eV 

647T 3 p 



1 - 



3^ 
p 2 



Sp 



P 2 \ 








(p - \i)pn 



(p + M J 



(41) 



4 

e p 



96tt 3 
eV 2 
32n 3 p 



1 + —T- log 1 

p z \ LM Z 

1 + 2^ + lo ( (P~ V)PA 

3p \(p + M ) 



(42) 



As in section ( [III B|) , the logarithmic divergence in q* corresponds to processes where the 
intermediate fermion state is soft in the second diagram of Fig. 8. In these processes, the 
incoming photon (boson) transfers almost all its energy to the outgoing electron. Since the 
outgoing electron is always above the Fermi level, this situation can occur only if p > fi (We 
assume for simplicity that \p — fi\ ^> eii). Thus the logarithmic divergence is present only for 
p > fi, as can be seen in eqs. (|4"1| ) and (|4"2"D. Then, both soft and hard momentum transfers 
contribute. 

The soft sector contribution corresponds to the emission of a soft antifermion (or to the 
absorption of a soft fermion) by a hard boson (Fig. 9, right). Since k! ~ p, the corresponding 
matrix element is approximately given by eq. (|20D for a Yukawa boson, and by 
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M± = eu(p, A) e T ■ ~fu(±q, A') 



(43) 



for a QED photon, where is the polarization vector. Since is transverse with respect 
to the photon momentum p, The operator et • 7 anticommutes with the Dirac operator 
P7° — p-~f and with the chirality operator 7 s . Thus, it simply changes u(p, A) into u(p, —A), 
up to a phase. The matrix element is the same as for the Yukawa interaction. The only 
difference is that the chirality A is conserved (A' = A), while it changes in the Yukawa theory 
(A' = -A). 

In both cases, the resulting soft sector contribution takes therefore the same form as the 
soft contribution for positron annihilation in eq. ( f23|) with an additional factor 2 for the 
final electron spins, and with an additional restriction on phase space from Fermi statistics, 
p — u>y. Since u is of order e/i, this restriction can be ignored as soon as the boson energy 
is not too close to the Fermi energy, i.e. if p — \x ^> e/i. In this condition, the damping rate 
is given by eq. (p?D , multiplied by a factor of 2. Adding up the hard and soft contributions, 
we find for the hard Yukawa boson 

r(p < A 



and for a hard photon 
V{p < 



4 

e p 


[l ^ 


' 1927T 3 






3p 


QAtt 3 p 


4 

e p 
96tt 3 




e 4 /i 2 " 


2>p 


32tt 3 p 



<•('-% . 

V mj(p + p.) t 



(44) 



+ log 



f 2py(p - p) 



(45) 



\ m 2 f (p + fi) 

Once again the cutoff g* has canceled out upon addition of the hard and soft contribu- 
tions. These results are valid only far from the Fermi surface, i.e. for \p — fi\ ^> e/i. The 
corresponding self-energy diagrams are displayed in Fig. 10. 



E. Gluon 

Three tree diagrams contribute to Compton scattering of a gluon. They are displayed in 
Fig. 11. The first, as shown below, yields a contribution of order e 3 /z; the contribution from 
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the other ones is subleading and is of the same order in e as in QED, i.e. of order e 4 /i or 
e 4 /zlog(l/e) (see previous subsection). 

The scattering process in Fig. 11, left, can be viewed as the emission of a soft gluon by a 
hard one (see Fig. 12). The matrix element is evaluated as follows. The three gluon vertex, 
coupling a hard gluon of color index a, momentum p and polarization e p to a hard and a soft 
gluon of color indices b and c, momenta p' and q and polarizations e p ' and e q respectively, is 

M = ef abc (e p ■ e p , (P + P') ■ e q + e p , ■ e q (-P 1 + Q) ■ e p + e q ■ e p (-Q - P) ■ e p >) , (46) 

where f abc is the SU(3) structure constant. The hard gluons are on-shell transverse gluons, 
whose polarization vectors satisfy e p ■ P = e p > ■ P' = 0. Therefore, in the limit q <C p, the 
last two terms in eq. ([R|) vanish. The remaining term is 

M = 2ef abc (e p -e p ,)(P-e q ). (47) 

This can be written in the form M = J p e^ where J p = 2ef abc (e p ■ ^ P ')P^ is the matrix 
element of the color current between the initial and final hard gluon states. In this form, it 
is analogous to the matrix element obtained in the case of the emission of a soft photon by 
a hard electron (see section [III C|) . 

The gluon damping rate is therefore given by an equation similar to eq. (|29|), with two 
minor modifications: the condition u < p — fi does not apply for a final gluon state, and the 
result must be multiplied by a factor N c /2 = fabcfabc/16 coming from the color degrees of 
freedom: iV c for the possible ways of transferring color to the soft gluons, and 1/2 for the 
symmetry factor. For N c = 3, the numerical value is 

T{p) = ^ 0.057 e 2 q D , 

= 0.019 JW f eV (48) 



IV. DISCUSSION 



Damping rates of hard particles in a cold ultrarelativistic fermion gas are at least of 
order e 3 higher than their energy: hence, one-particle excitations are narrow quasiparticle 

21 



states in the perturbative regime e < 1. The damping rates of the various one-particle 
excitations are displayed as a function of their momentum p in Figs. 14, 15 and 16 for the 
Yukawa interaction, QED and QCD respectively. 

We distinguish two categories of damping processes: for charged particles in gauge the- 
ories, the scattering process is essentially forward (8 ~ e) and the resulting damping rate 
is of order e 3 /x; for other particles, large angle scattering (8 ~ 1) contributes at least as 
much as small angle scattering, and the damping rate is of order e 4 /z or e 4 /zlog(l/e). Notice 
that the results derived in sections [III Q , |l fD] and |III E| are gauge invariant: all the matrix 
elements used in either the hard or the soft sector contributions are averaged over physical 
polarization states. 



A. Damping rates of charged particles in gauge theories 

Damping rates of charged particles are dominated by collisions with soft momentum 
transfer, for which medium effects must be taken into account: scattering takes place through 
the coupling of the elementary particle current with coherent plasma oscillations of the charge 
and current densities. The underlying classical structure is clear: the transition rate depends 
on the hard particle only through the associated current. We are in a situation where the 
hard particle motion is only slightly perturbed by the soft one: the gauge field behaves 
essentially as a classical field which couples to the current of the hard particles. 

The electron and quark damping rates are very similar, in the sense that they differ only 
by trivial color factors. On the other hand, the gluon and photon damping processes are 
essentially different: the dominant contribution to gluon damping involves the three-gluon 
vertex (see Fig. 11), and the damping rate T is of order e 3 /x, as for electrons and quarks. 
Gluon damping is therefore specifically non-abelian. The photon damping rate is smaller in 
magnitude, of order e 4 /ilog(l/e) or e 4 /z. 

It is interesting to note that the damping rates of charged particles are independent 
of the particle momentum. The only exception comes from electrons or quarks very close 
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to the Fermi surface, within an interval e/x from the Fermi level, where the damping rate 
decreases with the excitation energy \p — jj\ (section |111 C|) . Then, the longitudinal part 
of the interaction is screened at low momenta and leads to a width quadratic in \p — y\. 
The dominant term, shown in eq. fl3~T|), is the contribution from the transverse piece of the 
interaction, which is not screened in the static limit, and the width is linear in \p — 
This correlation between the range of interaction and the electron damping rate close to the 
Fermi surface is already well-known in the non-relativistic electron gas ( and In 
particular, our result of eq. (|37D agrees with the energy dependence of the imaginary part of 
the electron self-energy derived in ||20|| , if we set the Fermi velocity vp = pp/m to Vp — 1. 

Note that the calculation presented here is valid only for hard momentum excitations. 
The damping of soft, charged, excitations in gauge theories has very different properties: 
first, it is momentum dependent; second, the hard contribution is no longer negligible, but 
becomes of the same order of magnitude as the soft contribution, yielding a damping rate 
of order e 3 /ilog(l/e), instead of e 3 /z for hard excitations [IT) 



B. Damping rates of other particles 

Other particles include neutral particles (photons, Yukawa scalars) and fermions with 
a Yukawa coupling, for which the interaction <fnpi/} is not related to any conserved charge. 
For these particles, the hard contribution to the damping rate is at least of the same order 
of magnitude as the soft contribution. While the soft contribution involves medium effects, 
i.e. coherent effects, the hard contribution simply results from incoherent collisions: a 
quasiparticle excitation dies off by kicking electrons out of the Fermi sea randomly. The 
resulting damping rates are smaller in magnitude, of order e 4 // or e 4 /xlog(l/e). They are 
strongly momentum dependent, as can be seen in Figs. 14 and 15. 

The log(l/e) term comes from processes in which a massless fermion is exchanged, and 
the fermion propagator must be corrected for medium effects. These processes are specific 
to relativistic plasmas. As a result of the medium effects, the photon and the scalar boson 
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damping rates are strongly momentum dependent: they rise steeply near the Fermi energy 
(see Figs. 14 and 15), and are of order e 4 /i for p < p, and e 4 /ilog(l/e) for p > p. Once again, 
the soft momentum transfer process is almost classical in nature. Here, it is the fermionic 
soft field which acts as a classical field (recall that the dispersion relation of soft fermions is 
the same for both interactions) in which the hard particles move J2T|. 

Near the Fermi surface, the fermion damping rate in a Yukawa theory decreases quadrat- 
ically with \p — p\, in contrast with the electron and quark damping rates. This is a con- 
sequence of the fact that hard momentum transfers dominate and that in this sense, the 
interaction is short-ranged. 

Damping rates of soft excitations are of order e 3 p or e 3 /ilog(l/e), i.e. one power in e 



smaller than for hard momenta. The additional factor 1/e comes from kinematics [17 



C. Comparison with the high temperature case 

As noted in the introduction, ultrarelativistic plasmas have the same screening properties 
in the high density (T = 0) and high temperature (p = 0) limits. For neutral particles, we 
have seen that damping rates are of order e 4 /i or e 4 /i log(l/e) at T = 0. These damping rates 
are generally of order e 4 Tlog(l/e) at high temperature ||22j| . The log(l/e) factor comes from 



processes in which a soft fermion is exchanged, i.e. Compton scattering and pair annihilation 
at low angles. At T = 0, these processes are not always possible (see sections |III A| and QID) 



and the log(l/e) then disappears. Apart from this difference, the orders of magnitude of 
damping rates are the same in the hot and cold plasmas. 

For charged particles, the situation is very different. Naive perturbation theory yields 
an infrared divergent damping rate in the T = (p > 0) and p = (T > 0) limits. 
However, the level of divergence is different in these two cases, because screening at the 
one-loop level gives a finite damping rate, of order e 3 p, if T = 0, whereas in the high 
temperature case the damping rate is logarithmically divergent and of order e 2 T. In both 
cases, the damping processes are elastic collisions with the charges in the plasma. The 
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essential difference lies in Pauli blocking: at zero temperature, processes in which a boson 
with energy uj ~ e\x is exchanged have a phase space proportional to uj (only electrons very 
close to the Fermi surface participate in the collisions), whereas at high temperature, phase 
space is proportional to T . This factor uj both contributes a factor e and kills the divergence 
at uj = 0. 
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APPENDIX A: SPECTRAL DENSITIES OF SOFT MODES 

In this appendix, we recall how boson and fermion propagators at low momenta are 
modified by medium effects in an ultrarelativistic plasma. 



1. Soft gauge field 

In the Coulomb gauge, rotational invariance allows one to decompose the photon prop- 
agator into a longitudinal (L) and a transverse (T) piece [24] : 



D^(q, z) = D L {q, zy L el + D T {q, z) £ e»(q, A)e^(q, A) (Al) 

A=l,2 

where = 5$ and €^(q, A), A = 1,2, are spacelike unit vectors mutually orthogonal and 
transverse to q, therefore satisfying J2x=i,2 6 t(9> ^) e r(9' ^) = — lilj/Q 2 - F° r the bare 
QED interaction, the decomposition of the propagator DQ U (q, z) according to eq. ( [All ) gives 
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D l(q,z) = q 2 ; D ^(q, z) = z 2 - q 2 . 



(A2) 



For soft momenta q ~ e/i, the propagator is modified by medium effects: 

D uv(q, z ) = D oL(q, z ) - n^(g, z). 



(A3) 



To leading order in e 2 , the polarization tensor 11^ (q,z) is given by bubble diagrams cor- 
responding to the photon coupling to electron-hole intermediate states. Decomposing U^ u 



according to eq. (|Af|) , one gets (77 — > +0) |25| . |26 
f 



10 



IL L (q,w + i7i) = -l 

qD z q 



1 . u 

-5-n T (g,o; + 17]) = — + 



log 

1 u(q 2 -u 2 ) 



lo + q 



log 



w + q 



lu — q 



27T 



(A4) 
(A5) 



q D 2q 2 Aq 3 

Notice that the full one-loop self-energy contains also diagrams corresponding to electron- 
positron intermediate states, however their self-energies are a power of e smaller than those 
of eqs. (|A4j) , (|A5|) , sec |25| . [26|| . In eqs. (|A4|) , (|A5|) , qr> is the Debye screening momentum, 
given by the first entry of the following table. The resummed propagator obtained from 
eqs. (|A3|- |A5|) is drawn as a photon propagator with a "blob" (see for instance Fig. 2). 
The gluon propagator and gluon self-energy are diagonal in color indices, which we omit 



for brevity. In the strict Coulomb gauge, in the sense defined in the second reference of [] I 
the decomposition ( |A1| ) holds for the gluon propagator, and eqs. 



through ( |A5| ) remain 

valid. The Debye screening momentum qjj is given by the second entry of the following 
table. 



Theory 


qD 


QED 


ep/n 


QCD 


en^N~ f /{nV2) 



The spectral density p^ v (q,oj) = — 2 ImD'"'(q, u + irj), (7/ — > +0) can be decom- 
posed, like the propagator in eq. (|Al|), into longitudinal and transverse pieces, pL,r(q,^) = 
-2101^^(5,0; + it]). Their expressions are easily obtained from eqs. ( |A2| - [A5D . For free 
fields, they reduce to p L(q,u) = 0, p t(<?)^) = 27r<5(u; 2 — q 2 ){6{uj) — 6{—uj)): the only peaks 
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are at uj — ±q and correspond to transverse photons (gluons). For soft momenta, pi and 
Pt are modified by medium effects. The peaks of pr are shifted towards higher values of 
|a>|, and a peak appears in pl at the values of uj given by uj) — uj) = 0, corre- 

sponding to plasmon modes. Note that the real parts of Dj} T (q,uj) are even in uj for fixed 
q, so that the peaks always appear in pairs of opposite sign ±uj. In addition to these peaks, 
the spectral densities have a continuous part for |a;| < q coming from the imaginary part of 
the polarization in eqs. (|A4| -|A5|) , which corresponds physically to Landau damping: waves 
with |a;| < q lose their energy by accelerating fermions. 



2. Soft fermion 

The fermion propagator can be decomposed on a basis of spinors in the following way 
(once again, we omit trivial color indices for the quark propagator): 

G(p,z) = G+(p,z) ]T u(p,X)u(p,X) +G-(p,z) ]T u(-p,\)u(-p,X). (A6) 

A=-l,l A=-l,l 

In this decomposition, u(p, A) denotes a solution of the free massless Dirac equation {py°—p- 
~f)u(p, A) = with chirality A, normalized according to the relation J2\=-i,i u (p, A)w(p, A) = 
Pl°~P'l- Note that u(p, A) is a positive energy solution, while u(—p, A) is the corresponding 
negative energy solution with the same momentum p: it corresponds to a positron (or an 
antiquark) with momentum —p. An explicit expression of u(p, A) is most easily obtained in 
the chiral representation of Dirac matrices: 

u(p,+l) = ^>^ ) ^ u(p,-l) = ^^°J (A7) 

where p = p/p and (ftp is a two component spinor pointing in the direction of p, i.e. satisfying 
°" " P4>p — 4>pj normalized to unity (ftp<ftp = 1. Note that <ftp<ftp = (1 + cr ■ p)/2 

For the free Dirac propagator Gq 1 (p, z) = z^ — p ■ 7, the decomposition (|A6|) gives 

Gul(p,z) = 2p(zTp). (A8) 
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As expected, the poles of G + and G- are respectively the positive and negative energy 
solution of the free Dirac equation. For soft momenta p ~ ep, the fermion propagator is 
corrected by medium effects: 



G- 1 (p,z) = G^ 1 (p,z)-i:(p,z). 



(A9) 



The mass operator E(p, z) can be decomposed according to eq. (|A6|) , so that G^} = G Q ± 
E±. To leading order in e 2 , the components E± are given by (77 — > +0) | 



2m 2 



E + (p,w + i?7) = 1 



uj — p 
2p 



log 



u; +p 



uj — p 



— iix6(p 2 — u 2 



E_ (p, uj + irj) = —Re E + (p, —uj +irf) +i Im E + (p, —uj + irj) . 



(A10) 



The self-energy functions E-t only include the coupling of the soft external fermion to 
intermediate fermion states through the absorption or the emission of a boson. These terms 
are the dominant ones for a soft external fermion |^],[^|. The quantity m/ in eq. ( A10 ) is 
the quasiparticle rest energy and is given in the following table: 



Theory 


m f 


Yukawa 


e/i/(47r) 


QED 


e y u/(7rv / 8) 


QCD 


e/i/ (7r\/6) 



The spectral density p F (p,uj) = —2lmG(p,uj + irj) can be decomposed like the 
propagator in eq. (|A6|) , with G±(p,z) replaced by the corresponding spectral density 
p±(p,uj) = — 2 lmG±(p, uj + irj). For a free Dirac field, the spectral density reduces to 
p±(p, uj) = 2tx 5 {uj 2 —p 2 )6(±uj). Medium effects modify the spectral density for soft momenta. 
The position of the peaks of p + are the solutions of G^(p, uj) = 0. Using eqs. (|A8| - |A10|) , 
one obtains the dispersion relation in terms of the parameter x = uj/p: 

- 2 1 1, fx + 1 



P . 

m 2 x — 1 2 
uj = px. 



X 



(All) 
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There are two peaks for a given p. One with x > 1, which corresponds to the bare fermion 
state slightly shifted by its interaction with the medium. Furthermore, a second peak appears 
for x < —1. It corresponds to a new fermionic excitation called "plasmino" which has no 
counterpart in non-relativistic plasmas. Finally, p + has a continuous part in the region 
\uo\ < p, which corresponds to the fermionic analogue of Landau damping. The density p_ 
has the same properties, with uj replaced by — uj. 

APPENDIX B: SELF-ENERGY DIAGRAMS AND CUTTING RULES 

We show explicitly that the kinetic theory approach used in this paper is equivalent to 
the approach using the formalism of field theory. In field theory, the damping rate T(p) 
is defined from the imaginary part of the self energy. In the case of a fermion, using the 
notations of Appendix A, this relation reads (77 — > +0) p| 

r(p) = -2tr \jhnE(p,p + i7i)] /4p = -2 Im + irj). (Bl) 

Cutting rules (see for example |3(|) allow to express the imaginary part of a generic self- 
energy diagram as the rate of a scattering process, thus providing the equivalence with 
the kinetic theory approach used in this paper. We show in this appendix that the hard 
contribution to the damping rate corresponds to a two-loop self energy diagram where 
the loop momenta are hard (see Fig. 1) while the soft contribution corresponds to a one- 
loop self-energy diagram with a soft internal momentum (Fig. 2). We take the example 
of electron-electron scattering in QED. Our arguments can be easily extended to other 
scattering processes. 

We proceed as follows: we first show that the imaginary part of the one-loop resummed 
diagram depicted in Fig. 17 corresponds to the probability to emit a soft photon (second 
equality in Fig. 2). Then we show the equality displayed in Fig. 1, i.e. that the rate of 
electron-electron scattering corresponds to the imaginary part of a two-loop self-energy 
diagram. Note, however, that the interference term between the two Feynman diagrams of 
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electron-electron scattering (Fig. 3, right) is not included here. It corresponds to the imag- 
inary part of another two-loop diagram (Fig. 4, right). Finally, we show that the imaginary 
part of the one-loop resummed diagram corresponds to electron-electron scattering with a 
resummed interaction (first equality in Fig. 2). 



1. One soft loop 



The contribution of the diagram in Fig. 17 to the fermion self-energy is given by 

/da r+ioo A 7 
72^ Loo ^(-^GoiP-Q^p-^H^D^z,). (B2) 



In this expression, z p = fi + ix with x real, Go is the fermion propagator, which coincides 
with the free propagator ( |A8| ) for a hard fermion, and D is the soft photon propagator given 
by eqs. ([A3|-rA~5|) . We write the internal propagators using the spectral representations: 

50 dug p^(q,u q ) 



-oo 2,71 Zg — UJq 

+°° dujpi Pf{p'i Up 

x 2"7T Z n i — UJ r , 



with z q = ix, z p i = fi + ix and p' = p — q. A straightforward contour integration gives for 
the integral over z q : 

'+*°° dz q 1 1 1 + n{ui q ) - f(upf) 



ioo 2i7T Zq — UJq Zp — Zq — UJpi Zp — UJ q — UJpi 



(B4) 



where we have introduced the Bose-Einstein and Fermi-Dirac distribution functions n(u) 
and f{uj) which, in the limit T = 0, reduce to 

fM= e (»-MT + 1 = e b-<*)- ( B5 ) 

After analytic continuation of z p to uj p + it], the imaginary part of eq. ([B4]) becomes 

7l5(uJp -UJg- UJ p >) (1 + n{Uq) - f{uJ p ')) . (B6) 
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The imaginary part of the self-energy can thus be obtained from eq. ( |B2| ) through replacing 
the internal propagators by their spectral functions and the Matsubara frequencies z by real 
frequencies u, and multiplying by the occupation factor from eq. (|F36|): 

-21mV(p,p + ir,) = J ^(er)PF(P-Q)(eY)p,u(Q) 

x (1 + n(u q ) - f(u p - u q )) , (B7) 

where we have introduced the four vectors P = (p,p) and Q = (u q , q), and u p = p. 
Note that the occupation factors can be rewritten as 

1 + nK) - /(ay) = (1 + n{u q )) (1 - /(ay)) + n{u q )f{u p ,). (B8) 

The two terms correspond to the amplitudes for the direct and inverse process in Fig. 2 
(right). If Up > ii (particle excitation), only the direct process contributes, when < uj q < 
u p — \l. If Up < \i (hole excitation), on the other hand, only the inverse process contributes, 
when \i — Up < u q < 0. Note that in this last case, the occupation factor is nf = — 1; 
however, the boson spectral function also has an opposite sign for u q < (see eqs. ( S3} - [^SD ) 
so that the global sign is unchanged. 



Decomposing the spectral functions according to eqs. ( |A1| ) and ("KB), and taking into 
account that p+(p, u) = 2n5(u 2 — p 2 )6{u) for a hard fermion (the spectral function is the 
same as in the vacuum, and p_ does not contribute since u p — u q > 0), one obtains from 

eq. my- 

-2u(p,\)ImZ(p,p + ir ] )u(p,\)= [ ^%2tt5((P - Qf) £ Pi (Q)£|M/ 

J j=L,T A' 

x(l+n(u q )- f(u p -u q )). (B9) 
where Mj is the matrix element of the transition process in Fig. 2 (right): 

Mj = u(p', A'X-ieTXp, X)e^(q). (BIO) 
Comparing eq. ( P9|) with eq. (0), one concludes 

r(p) = -2u(p,X)lmi:{p,p + ir])u(p,X)/2p, (Bll) 

= -2 tr [p 1 ImS(p, p + irj)] /4p. (B12) 
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2. Two hard loops 



We now turn to the second step: we show that the imaginary part of the two-loop 
diagram on the right of Fig. 1 corresponds to the rate of electron-electron scattering with 
the diagram on the left of Fig. 1. The contribution of the two-loop diagram to the self- 



J-ioo ^\G*{k,z k ){e^)G,{k + q ,z k + z q ){eY)] (B13) 



energy is given by a formula analogous to eq. (|B2|), where D is replaced by a photon line 
with a fermion loop insertion, i.e. by D UD , with 

a r dzt ' 

I fl— too 

with z q on the imaginary axis. Therefore, the imaginary part of the two-loop diagram is 
given by eq. (|B~7|), in which the boson spectral function p(q,u q ) = —2ImD(q,ui q + irj) is 
replaced by —2 lm(D UD ). Now, the imaginary part of D vanishes because both fermions 
are on mass shell {\uj p \ = p and \u p — u q \ = \p — q\), which implies \u\ < q. Thus we can 
write lm(D UDQ) = D (lmU)D . To calculate Imll, we follow the same steps as for ImE 
in eq. ( p2|) . Using the spectral representation to write the internal propagators, the integral 
over Zk can be calculated easily {k' = k + q): 

r +io ° dzk_L 1 = I'M - I'M) 

Jfi-ioo 2ilT Zk — LOk Zk + Z q — LUk< Z q + UJk — U)}.' 

After analytic continuation of z q to u q + irj, the imaginary part of this equation becomes 

- 1l5(uj q + UJk~ UJ k ') (I(uJk) - fM)) ■ ( B15 ) 

One thus obtains 

(2^)3 ^f tr [Pi?(fc ' Uk) ( e ^ PF ( k + w * + u i)M)} 
x(fM-fM))- (B16) 



Replacing p in eq. (|B7|) by — 2D (lmH)D , using eq. (|B16| ), and decomposing the fermion 
spectral functions according to eq. (|A6|), one obtains 



r(p) = i / S3i f ^4 M)5((p - Qf) M5M) Mm + Qf) \m\ 2 



2pJ (27r) 4 J (2tt) 4 

x (f(u k ) - fM))(l + n(u q ) - fM) • (B17) 
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where K = (uk, k). Notice that, from K 2 = 0, the only possible solution for u>k is oof. = k, 
only electrons are initially present in the Fermi sea. The matrix element M is given by 

M = Ea',«a' J p',\>;p,\ J k',K>;k,K,iO-/Q 2 ) with Jp, tX <;p,\ = eu(p' , X')Yu(p, A) and corresponds 
indeed to the direct contribution to electron-electron scattering (Fig. 1, left). Introducing 
the four vectors K' = K + Q and P' = P — Q, integrating over Q and making use of 
the identity / \d A K) / (2ii) A 27i5(K 2 )Q(K ) = J drk, one finds eq. (f|), up to the statistical 
factors. Thus, we only need to check that the product of the phase space factors eq. ( |B6|) 
and ( p3 1 5|) corresponds to electron-electron scattering. For this purpose, we note that the 
energy conservation u q + Uk = ujk> implies the following relations between the statistical 
factors: 

™K) (fM - f{uk>)) = (i - fM) /(wfc/) 

(1 + n{u q )) (fM - f(u k ,)) = (1 - fM) f(u k ) (B18) 

Using these equations together with eq. (P8|), one obtains immediately the phase space factor 
under the form 

(1 - fM))fM){\ - fM)) + /(uy)(l - fM)fM)- (B19) 
The two terms correspond to the amplitudes of the direct and inverse process, as expected. 

3. Screened interaction 

Finally, we show that the imaginary part of the one loop resummed diagram corre- 
sponds to the probability of electron-electron scattering with a screened interaction (first 
equality in Fig. 2). This is a straightforward extension of the previous result. We 
start from eq. (|B7|). The spectral function of the resummed photon line, p, is given by 
p = — 2IulD = 2\D\ 2 (ImD~ 1 ). Using eq. (|A3|) and the fact that IulDq 1 = 0, as discussed 
above, one obtains p = — 2\D\ 2 (lmU). Thus the only difference with the previous case is that 
the free photon propagator Dq is replaced by the resummed photon propagator D, which 
includes the screening effects. 
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FIGURE CAPTIONS 



Figure 1: Left: Tree diagram for e~e~ scattering. Right: Two-loop self-energy diagram. 
The imaginary part, obtained by cutting the diagram through the fermion loop, corresponds 
to the amplitude on the left, squared and integrated over phase space. 

Figure 2: Left: Diagram for e~e~ scattering with a screened interaction. The resummed 
photon propagator is indicated by a blob. Middle: One-loop self-energy diagram with a 
resummed photon propagator. Right: Emission of a virtual soft photon by a hard electron, 
corresponding to small angle e~e~ scattering. 

Figure 3: Tree diagrams for M0ller scattering. 

Figure 4: Two-loop self-energy diagrams corresponding to the scattering processes depicted 
in Fig. 3. The diagram on the left gives the direct and exchange contribution while the 
diagram on the right is the interference term. 

Figure 5: Tree diagrams for (a) Bhabha scattering and (b) pair annihilation. 

Figure 6: left: Tree diagram for e + e~ annihilation. As in Fig. 2 (left), the exchange diagram 
is negligible when the momentum carried by the internal propagator is soft. Right: Emission 
of a soft virtual positron (or absorption of a soft electron) by a hard positron, representing 
the contribution of soft momentum transfers to the process on the left. 

Figure 7: Self-energy diagrams corresponding to positron scattering and annihilation. The 
three two-loop diagrams correspond to the hard contribution: the first diagram corresponds 
to direct and exhange terms in Bhabha scattering (Fig. 5a), the second to direct and exchange 
terms in pair annihilation (Fig. 5b), while the third diagram, which can be cut in two different 
ways, gives the interference terms of both processes; finally, the one-loop diagram on the 
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right corresponds to the soft contribution to pair annihilation. 

Figure 8: Tree diagrams for electromagnetic and Yukawa Compton scattering. 

Figure 9: Left: Leading diagram for Compton scattering in the limit of soft exchanged 
momenta. Right: Emission of a soft virtual positron (or absorption of a soft electron) by a 
hard photon, corresponding to the process on the left. 

Figure 10: Self-energy diagrams corresponding to Compton scattering. The three two-loop 
diagrams correspond to the hard contribution: the first and second diagrams correspond 
respectively to the first and second processes in Fig. 8, while the third diagram is the inter- 
ference term; finally, the one-loop diagram on the right corresponds to the soft contribution. 

Figure 11: Tree diagrams of gluon Compton scattering. The first of the three diagrams 
dominates when the momentum carried by the internal gluon is soft. 

Figure 12: Emission of a soft virtual gluon by a hard gluon, representing the contribution 
of soft momentum transfers to the process depicted in Fig. 11. 

Figure 13: One-loop self-energy diagram contributing to the gluon damping rate through 
Compton scattering at small angles. 

Figure 14: Damping rate of hard one-particle excitations in the Yukawa theory as a function 
of their momentum p. We have chosen for e the numerical value e 2 /47r = 1/137. Full line: 
fermion (p > fi, eq. (|T2"D) and hole {p < fi, eq. (0)) excitations; long dashes: antifermion 
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(eqs. ( fP7|) and (|28|)); short dashes: boson (eq. ([44])). 
Figure 15: Same as Figure 14 for QED interaction. 

Figure 16: Same as Figure 14 for QCD interaction, with Nf = 2 flavors. The numerical 
value of the coupling constant is the same as in Figs. 14 and 15. 

Figure 17: One loop resummed diagram contributing to the electron or quark damping rate. 
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